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1. Introduction 



The goal of this paper is to construct a QCD motivated model for the strong interactions at high energy, 
and to elucidate the model's predictions and implications at LHC and Cosmic Rays energies. 

The difficulties and challenges of a theoretical approach to the strong interaction are well known: a 
qualitative understanding of the confinement of quarks and gluons in QCD. For high energy scattering the 
situation is even more difficult than for the hadron structure, since approximate methods such as QCD 
sum rules and/or effective theories, as well as the lattice QCD approach, cannot be used to calculate the 
high energy amplitudes. Today, and for the past four decades, the accepted method of describing soft 
interactions at high energy, is the phenomenology based on the soft Pomeron and secondary Reggeons (see 
Refs. [1-3] for details). All parameters related to the Pomeron and Reggeons, which are assumed to be 
simple poles in the J-plane, such as the intercepts and slopes of the trajectories, the vertices and their 
dependence on the impact parameter, have to be deduced from the experimental data. However, the key 
problem of all approaches based on the soft Pomeron hypothesis, is that there is no theory which can 
specify what kind of interactions between Pomerons have to be taken into account, as well as the values 
of the multi-Pomeron vertices. Due to this problem we are doomed to build basically unreliable models, 
since the only criteria is that they describe the experimental data (see Refs. [4-7]). 

In this paper we attempt to overcome this difficulty by using high energy perturbative QCD approach. 
At first sight, this appears unrealistic, since the pQCD approach is based on the smallness of the running 
QCD coupling at short distances, while the high energy interaction is a typical example of long distance 
non-perturbative QCD. We wish to question this widely held prejudice. Indeed, the only microscopic expla- 
nation for the Pomeron structure is given in the partonic approach [8] , in which the slope of the Pomeron 
trajectory is related to the mean transverse momentum of the exchanged partons {a'p oc 1/ < pt >, where 
< pt > is the mean parton momentum). The commonly held view in high energy phenomenology, is that 
a'p = 0.25 GeV~ 2 [9], which implies that in the Pomeron the typical momenta of partons are high (espe- 
cially if you compare this value with the value of a'p, = 1 GeV~ 2 for the secondary Reggeons). In the next 
section, we extend our output up to the Planck mass, and will show that the fit to the experimental data 
based on our model with power-like dependence of the Pomeron-proton vertices on energy, validates the 
consistency of our calculations with unitarity. The momentum transferred t behaviour leads to a'p < 0.02 
GeV~ 2 . This result, together with the fit of Ref. [7], where the data were fitted with a'p = 0, lends support 
to our assumption that, the typical parton momentum is large (approximately < pt >= 1/y^a'p > 7 GeV). 

Therefore, the running QCD coupling as = 7r/61n ^< p 2 > / h^cD*j ^ 1 (approximately 0.18), and we can 
consider it as our small parameter, when applying perturbative QCD estimates to the Pomeron-Pomeron 
interaction vertices. 

The theoretical scheme that emerges from such an approach will be discussed in section 3. Using per- 
turbative QCD we select the essential Pomeron-Pomeron vertices, and develop the method for summation 
of all these diagrams. We also give the partonic interpretation of our approach, and develop an approxi- 
mate method for the solution of the long standing problem of the high energy asymptotic behaviour of the 
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scattering amplitude. Based on this analysis we derive formulae for all cross sections taking into account 
the Pomeron interactions. 

In section 4 we fit the available experimental data using our formalism. This fit is based on an updated 
data base which includes the published p-p and p-p data points of atot, the integrated values of cr e ;, a s d, 
add and the forward elastic slope B e [, in the ISR-Tevatron energy range, p = f^a'^t^o^) anc ^ ^he forward 
slope of the SD and DD final states are predictions of the model. The additional data that we have used 
to validate our present model, are the t-dependence of the differential elastic cross section, and the mass 
dependence of the single diffractive production. We successfully describe this data. Based on this fit we 
give reliable predictions of the quantities to be measured at the LHC cm. energy of lATeV. Our output 
also covers the broad Cosmic Rays energy range up to the GZK limit. 

The fifth section is devoted entirely to a calculation of the survival probability ((| S 2 |)) for exclusive 
central diffractive Higgs production at the LHC, and a discussion on the reliability of these calculations. 
We confirm the tendency in which (| S 2 |) becomes small (less that 1%), noticed in our previous paper [6]. 
From a practical point of view, this is the most salient feature of this paper. 

In the conclusions, we compare critically our results with other approaches, and summarize our findings. 
We list a few experimental signatures, which should enable us to differentiate between alternative theoretical 
options and phenomenological models. 

2. The two channel model and the value of ct'^, 
2.1 GLM two channel model 

The GLM two channel model has been described in our previous publications (see Refs. [4,6,10-12], and 
references therein). In this formalism, diffractively produced hadrons at a given vertex are considered as 
a single hadronic state described by the wave function typ, which is orthonormal to the wave function Vl/^ 
of the incoming hadron (proton in the case of interest), < Vf/il^ d >= 0. We introduce two wave functions 
tpi an d ip2 that diagonalize the 2x2 interaction matrix T 

A,k =< tpi tpk\T\ipi' ipk' >= A i,k Si,i> Sk,k>- (2.1) 
In this representation the observed states are written in the form 

^ = aVi + ^2, (2.2) 

^D = -/3tpi + atp 2 , (2.3) 
where, a 2 + fi 2 = 1. Using Eq. ([2,1|), we can rewrite the unitarity constraints in the form 

Im A i>k (s, b) = \A i>k (s, b) | 2 + Gf k {s, b), (2.4) 

where Gft is the contribution of all non diffractive inelastic processes, i.e. it is the summed probability for 
these final states to be produced in the scattering of particle i off particle k. 
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A simple solution to Eq. (|2.4| ) has the same structure as in a single channel formalism, 

A lik (s, b) = i ( 1 - exp (-*M?lR)) , (2.5) 



Gfl(s,b) = l-exp(-n itk (s,b)). (2.6) 



From Eq. (2.6) we deduce, that the probability that the initial projectiles (i,k) reach the final state 



interaction unchanged, regardless of the initial state rescatterings, is Pf k = exp (— Oj^(s, b)). 
For the opacities we use the expression 

fti,k{s,b) = gig k l—j S {b-mi,m k -a'jp\n(s/s )) (2.7) 

which differs from the expression that we used in our previous models. The profile function 

S (b, a'jp] mi, m k ; ]h(s/sq)) at s = so, corresponds to the power-like behaviour of the Pomeron-hadron 

vertices, 

S(b;m i ,m k ;,a' 1P ln(s/s ) = 0) = J 9i {q) g k (q) e^, (2.8) 

with a normalization J d 2 b S (6; rrii, m k ; , a'jp ln(s/so) = 0) =1. 

In this paper we choose 

(1 + q z /mf) z 

The arguments we list in favour of this choice are: (i) pQCD leads to gi(q) — ► ct'g(q)/q 4: [13]; (ii) in some 
models (for example in the constituent quark model) gi{q) is equal to the electro-magnetic form factor 
of the proton, which has the form of Eq. ( |2.9| ); and (iii) Eq. ( |2.9D reproduces the experimental elastic 
differential cross section for p-p (p-p) data in the range of t, up to 1 to 1.5 GeV 2 [9], while a Gaussian 
parametrization fits the data only for very small t < 0.1 GeV 2 . 

Using Eq. fl2,8|) we obtain for S (b; rrii, m k ; , a'p ln(s/so) = 0), 

» S (b;m,i,m k ; ,a']pln(s/so) = 0) = (2-10) 



(1 + q 2 /m 2 ) 2 (1 + q 2 /m 2 ) 2 

{Ami m k (K (mj6) - K (m k b)) + (m 2 - m 2 k )b {m k K x (m^) + m^i (m k b))} . 



4tt (mf — m 2 ,) 3 

For energies s > Sq, we need to take into account the observed shrinkage of the diffraction peak. To this 
end we replace gi(q) g k (q) by g%{q) g k (q) exp (— a'jp m(s / so)q 2 ) in Eq. ( pjjj) . To simplify our calculations we 
replace 

2 

m 2 => m 2( s ) = m ; (2.ii) 

1 + a^p ln(s/s )/4m 2 
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It is easy to check that gi(q; mi) <?&(<?; mk) exp (—a'p ln(s/so)q 2 ) and gi(q\ rrii(s)) <?&(<?; TUk(s)) have the same 
behaviour for a'p ln(s/so)<? 2 <C 1. When a'p ln(s/so)g 2 3> 1 these two expressions are different. Note 
that the Regge factor exp (— a'p ln(s/so)§ 2 ) cannot be justified in this kinematic region. In the region of 
a'p ln(s/so)q 2 ~ 1 it is preferable to use the Regge factor. We assume that a'p is small enough, so that 
this region gives a negligible contribution to all experimental observables. 

In general, we have to consider four possible re-scattering processes in Eq. fl2,l|) . However, in the case 
of p-p ( p-p) the two non-diagonal amplitudes are equal A\^ 2 = A 2j i, and we end up with three rescattering 
amplitudes. These amplitudes are presented in our two channel formalism in the following form [4, 11, 12] 

a e i(s, b) = i{a A A hl + 2a 2 (3 2 A^ 2 + (3 4 A 2i2 }, (2.12) 

a sd (s, b) = iap{-a 2 A hl + (a 2 - f3 2 )A h2 + p 2 A 2<2 }, (2.13) 

a d d = ia 2 p 2 {A hl - 2A 1>2 + A 2>2 }. (2.14) 

It should be stressed that in this approach diffraction dissociation, appears as an outcome of the Good 
and Walker mechanism [14] (G-W) or, in other words, the elastic, single diffraction and double diffraction 
processes occur due to elastic scatterings of ipi and tjj 2 , the correct degrees of freedom. 

The corresponding cross sections are given by 



<r to t00 =2 J d z ba el {s,b), (2.15) 

a el (s) = J d 2 b\a el (s,b)\ 2 , (2.16) 

a sd (s) = J d 2 b\a sd (s,b)\ 2 , (2.17) 

(s) = J d 2 b\a dd (s,b)\ 2 . (2.18) 



2.2 a'p ► 



Using Eq. (2.1) - Eq. ( |2.7D we fit the experimental data, so as to find the value of a'p. However, most of 



the data are available at rather low energies W = y/s = 20 to 600 GeV, where the contributions of the 



secondary Reggeon have to be included. For this we replace fij £ (s,b) given by Eq. (2/7) by the sum of 
Pomeron and Reggeon contributions, 

n i>k (s,b) = n? k (s,b) (seeEq. Q) + (2.19) 

with 

n-(sb) - 9 l^Ml e J O (2 20) 
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P 


a'jp 


9i 


92 


777-1 


m 2 


0.120 


0.46 


0.012 GeV' 2 


1.27 GeV- 1 


3.33 GeV' 1 


0.913 GeF 


0.98 GeF 


A K 


P 


a 'n 


9? 


9f 


^0,1 


X 2 /d.o.f. 


-0.438 


0.46 


0.60 GeV- 2 


4.0 GeV' 1 


118.4 GeV' 1 


4.0 GeV- 2 


0.87 



Table 1: Fitted parameters for two channel (eikonal) model 



The signature factor is 



corresponding to a Reggeon trajectory am(q) 



"L _|_ e «7rajR(g) 
sin (7ra K (<7)) 
1 + Ajr + a' R q 2 , and 



so 



(2.21) 



(2.22) 



gf'jB.Qi , Ajr and are fitted parameters. For the Regge sector we know the natural values for Ajr and 
a'jft from the behaviour of the Regge trajectory at q 2 < (t > 0) , since all hadrons lie on these trajectories, 
A K « -0.5 and 1 GeV- 2 . 

Eq. ( p. 19 ) - Eq. ( p. 22 ) specify our model. In this model we neglect the interactions between Pomerons, 
but include eikonal type rescatterings. This means that we assume only a G-W origin for all quasi-elastic 
processes. The fitted parameters are shown in Table 1. The quality of the fit is illustrated in Fig. [I] and 
Fig. H|, where our fit is shown by dashed lines. The fit is very good with x 2 /d.o.f. = 0.87. However, 
this parametrization fails to fully describe the cross sections for single and double diffractive production, 
yielding values for these cross sections which are approximately two times smaller than the experimental 
values. 

From this fit we have two conclusions: i) we cannot describe the diffractive production in the framework 
of the G-W mechanism if we assume an exponential parameterization for the proton profile function; and (ii) 
the value of a'jp turns out to be very small a'jp = 0.012 GeV~ 2 . As we have discussed in the introduction, we 
conclude from the smallness of a'jp, that the hard processes which occur at short distances are responsible 
for the Pomeron structure. 



3. Pomerons Interactions 
3.1 QCD input 

Our main hypothesis is that the typical distances for a Pomeron exchange are short, and we can use pQCD 
as a guide for building a theory for Pomeron interactions, based on the small value of a'jp, the Pomeron 
slope that we obtained in the previous section. We recall that the value of a'jp in the parton model [8] can 
be written as 

R 2 (s) = (b 2 ) n = (Ab 2 )n = -J^ p \n(s/s ) = a'jpHs/s ) (3.1) 

< Pt > 



- 6 - 



: o tot (s)(mb) 
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.... 1 .... 1 .... 1 .... 1 .... 1 ... . 


1 .... 1 .... 1 



2.5 3 3.5 4 4.5 5 5.5 

log 10 (s/s Q ) 



6.5 



Figure 1: Energy dependence of a tot- The solid line 
shows the fit with taking into account all Pomeron 
interactions while the dashed line corresponds to two 
channel (eikonal) model. 




4 4.5 5 5.5 

log 10 (s/s ) 



6.5 



Figure 2: Energy dependence of the slope for the 
differential elastic cross section. All notation are the 
same as in Fig. [I] 



where n is the number of partons at rapidity Y = ln(s/so) and p is theb density of partons in a unit of 
rapidity. In Eq. ( |3.1|) the radius of interaction, is determined as the average distance in impact parameter 
space that a 'wee' parton can reach at given enegy. 

In Eq. ( |3.ip R 2 (s) is the radius of interaction at energy W = \fs, which is the distance in impact 
parameter b space, that is reached by partons in a two dimensional diffusion, after n steps. (Ah 2 ) is the 
mean displacement during one diffusion step, which can be estimated using the uncertainty principle as 
(Ab 2 ) = 4/ < p 2 >. < pt > is the average parton transverse momentum which we wish to evaluate. The 
number of diffusion steps is equal to the number of partons, since at each step one parton is emitted. At a 
given energy the average number of emitted partons is equal to pm(s/so), where p is the parton density in 
units of rapidity. Comparing n with the hadron multiplicities, we conclude that p ~ 1 or larger. Therefore, 
Eq. (U) leads to < p\ > 20GeV for a'p = 0.01 GeV~ 2 . The typical value of the running QCD coupling, 
which corresponds to this value of the parton momentum, is as ~ 0.15, which is small enough to be used 
in pQCD estimates. 

In our procedure we only use general results of the perturbative QCD approach to high energy scattering 
[15]. Consequently: 

1.) In the leading order approximation of pQCD, only a Pomeron splitting into two Pomerons and 
two Pomeron merging into one Pomeron, should be taken into account [16,17], while all other vertices are 
small. Therefore, using this input from pQCD, we restrict ourselves by summing Pomeron diagrams with 
triple Pomeron vertices only. 



-7- 



2) Since Aafjpln(s) <C 1 over the entire kninematical range, we have investigated, we can neglect the 
value of a'jp, and consider the theory with a'jp = 0. 

3) We obtain the anticipated values for all ingredients of the Pomeron interaction approach: the 
intercept of the Pomeron Ajp above unity oc as, and the triple Pomeron vertex coupling g^jp oc a|. 

The theory that includes all the above ingredients can be formulated in functional integral form [17], 

Z[$,$+] = J D§ D$ + e s with 5 = S + Sj + S E , (3.2) 

where Sq describes free Pomerons, 5/ corresponds to their mutual interaction and Se relates to the inter- 
action with the external sources (target and projectile). Since a'jp = 0, Sq has the form 



So = J dY$+(Y) j-A + A j 



(3.3) 



Si includes only triple Pomeron interactions and has the form 

Si = 93P J dY {$(y)$+(y)$+(y) + h.c.} . (3.4) 
For Se we have local interactions both in rapidity and in impact parameter space, 

r 2 

S E = - I dY^ {*(Y) 9i (b) + $+(y) , (3.5) 

i=i 

where gi(b) stands for the interaction vertex with the hadrons at fixed b. 

In the next sections, after specifying this theory, we solve it. Indeed, this theory as any theory of 
Pomeron interactions, is written in such a way that an high energy amplitude satisfies t-channel unitarity. 
However, s-channel unitarity remains a problem. In the next subsection we will change the interaction 
term (Si), in such a way, that the theory will have a clear partonic interpretation, and will, also, satisfy 
s-channel unitarity. 

3.2 Generating function and a partonic interpretation 

To find a reformulation given by the functional of Eq. ([T^), we consider a system of partons* that can 
decay and merge: one parton to two partons, and two partons into one parton, with probabilities T(l — > 2) 
and r(2 — * 1), respectively. For such a system of partons, we can write a simple equation. Indeed, let P n {y) 
be the probability to find n-parton (dipoles) with rapidity y in the wave function of the fastest (parent) 
parton (dipole), moving with rapidity Y > y. For P n (y), we can easily write down a recurrence equation 
(see Refs. [20,21]) 

d P (v) 

_ = r(1 _> 2 ) {- n P n + (n-l)P n _i} + r(2->l) {- n (n - 1) P n + {n + l)nP n+l }. 

oy 

(3.6) 



*The partons for high energy QCD are the colourless dipoles, as was shown in Ref. [18]. 
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In each bracket the first term on the r.h.s., can be viewed as a probability of a dipole annihilation in the 
rapidity range (y to y — dy) (death term). The second is a probability to create one extra dipole (birth 
term). Note the negative sign in front of dP n (y)/dy. It appears due to our choice of the rapidity evolution 
which starts at the largest rapidity y = Y, of the fastest dipole and then decreases. The first two terms are 
responsible for the process of parton decay, while the last two terms describe the contribution of partons 
merging. 

It is useful to introduce the generating function [18, 19, 21] 

Z(y, u) = p n(v) u n (3.7) 

n 

At rapidity y = Y there is only one fastest parton (dipole), which is P\{y = Y) = 1, while P n >i(y = 
Y) = 0. This is the initial condition for the generating function 

Z{y = Y)=u. (3.8) 

At it = 1 

Z(y, u = 1) = 1, (3.9) 
which follows from the physical meaning of P n as a probability. 

Eq. (p.6|) can be rewritten as an equation in partial derivatives for the generating function Z(y,u), 



dZ(y,u) „, , , dZiy.u) „. , . d 2 Z(y, u) 

— j^—t = -r(i->2)t*(i - u) ' + r(2->i)it(i - u 1 ■ (3.io) 

ay ou a A u 



The description of the parton system given by Eq. ( 3.10 ), is equivalent to the path integral of Eq. (3.2). 
Indeed, the general solution of Eq. (|3 . 1 0|) has a form 



Z(Y;u) = e H{u)(Y - Yo) Z(Y ;u), (3.11) 

with the operator H defined as 

H (u) = - T(l - 2) u(l - u) — + T(2 - 1) u(l - u) jQ-y, (3.12) 

and 

Z(Y ;u) = e ^ h){u - 1] . (3.13) 
We introduce operators of creation (a + ) and annihilation (a) that satisfy [d, d + ] = 1, at fixed Y. 

d 

a = —,a + = u. (3.14) 
ou 

In this formalism 

H = -T(l -> 2)d + (1 -a + )a + r(2 -»■ 1) d+ (1 - d+) a 2 , (3.15) 
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and the initial state at Y = Yq is denned as 

|y > = e 9(&X« + -i)|o >, (3.16) 

with the vacuum defined as a|0 >= 0. 

The theory with the Hamiltonian of Eq. ( |3.15| ) , has an equivalent description using the path integral 
of Eq. (|3.2| ) (see the detailed derivation in Ref. [24]) with 



S= /( $ + H($ + + 1,-$)) dY = I dY ( — $ (3.17) 



dY v / J V dY 

- r(i -> 2)$+$ + r(i -► 2)$+($) 2 + r(2 -» i)(<i> + ) 2 $ - r(2 -► i)($+) 2 ($) 2 ) . 

Comparing Eq. (fTT7D with Eq. - Eq. (|3l|), one can see that they are similar if we put T(l — > 2) = A 
and r(2 1) = gj F /A. However, Eq. has an additional term (-T(2 -> l)($+) 2 ($) 2 ) which 



describes the two Pomeron to two Pomeron transition (four Pomeron interaction). This term ensures that 
our approach satisfies s-channel unitarity (see Ref. [25]), where it is shown that without this term there is 
no probabilistic interpretation of the Pomeron interaction theory, and s-channel unitarity is violated. We 
also need to renormalize $ + — * {gsip/ A) $ + and $ — * (A/ g%jp) <£. 

Using the functional Z, we find the scattering amplitude [22], using the following formula: 

N(Y) = ImA el (Y) = f]t±t dnZ< f> u) \ u=lln (Y = Y ,b), (3.18) 

n=i 

where ^ n iX = ^o,b) is the scattering amplitude of n-partons (dipoles) at low energy. These amplitudes 
depend on the impact parameters which are the same for all n partons, since a'jp = 0, and we neglect the 
diffusion of partons in impact parameter space. Eq. ( |3.18| ) corresponds to the partonic approach [8,23], 
in which a high energy scattering can be viewed as a two stage process. The first stage is a development 
of the partonic wave function, which we consider by introducing the generating function Z. In Eq. ( |3.18 ), 



the second stage is the interaction of the lowest energy partons ('wee' partons) with the target, which is 
described by the amplitudes 7 n (^ = Yq, b). Assuming that there are no correlations between the interacting 
partons (dipoles) at low energy, we can consider ^y n (Y = Yo,b) = "fi(Y = Yo,b) [22], which we include, 
choosing S# in the form of Eq. (|3.5| ) and of Eq. ( |3.16| ). 

The generating function approach given by Eq. ([D]), Eq. ( |3.10[) and Eq. flOf ), has the advantage that 



it can be solved analytically (see Ref. [26]). This solution leads to a constant cross section at high energy, 
while the interaction without the four Pomeron term, decreases at high energy [27]. This fact emphasizes 
the importance of s-channel unitarity, in finding the asymptotic behaviour of the scattering amplitude at 
high energy. Having an exact solution we are able to develop approximate methods, which we can check 
against the exact solution. 
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Figure 3: The different contribution to the first enhanced diagram and its renormalization procedure 



3.3 Improved Mueller-Patel-Salam-Iancu approximation 



Calculating the high energy amplitude in the generating function approach, we use the approximation 
which allows us to write a simple analytical formulae for the physical observable. The main idea of this 
approximation, which we call the improved Mueller-Patel-Salam-Iancu approximation, is the following: we 
claim that at high energy, in the kinematical region 



Y < 



A 2 
9zjp 



1 

J 

7 



(3.19) 



only large Pomeron loops, with a rapidity size of the order of Y, contribute to the high energy asymptotic 
behaviour of the scattering amplitudes. This approximation has been discussed in detail in Refs. [28-30], 
and here we illustrate this idea using the example of the first Pomeron loop diagram given in Fig. H Using 
Eq. (|3.2[) or the generating function approach, we obtain the contribution of this diagram in the form^ 



^In Eq. (3.20) we redefine the hadron-Pomeron vertices denoting them by g(b) -^y. We will clarify our reason for doing so 
below. For simplicity we use g(b) for this vertex, neglecting the different hadronic states in our two channel model. 
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A (Fig.® = (3.20) 

rvi 

/ dy 2 G(Y-y 1 )G 2 (y 1 -y 2 )G(y 2 -0) 
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l'AzpY _ _J_ e AjpY _ Y c Ajpy 
A|p Apom 



„A v omY~ 

C - T 7 ^1 ^ — r xj J. J B F 

A R Y .,2 JApYl 



— 9 2 {b){ lR e^ Y - ^e 2 ^ Y }, (3.21) 

with 7tj = 7 + 7 2 , and the renormalized intercept of the Pomeron Are = Ajp+T(2— > 1) = Ajp + Ajpj. 
In the kinematic region of Eq. ( p. 19 ) we can neglect the renormalization of the intercept, since in this region 



A;p7^ < 1. 

The general procedure for summing large Pomeron loops was suggested by Mueller, Patel, Salam and 
Iancu (MPSI) in Ref. [28]. In this approach, the scattering amplitude is calculated using the unitarity 
constraints in the t-channel (assuming that the amplitudes at high energy are purely imaginary, namely 
N = ImA), 

N([...]\Y) = N([...]\Y-Y';P^nP) N([. . .]\Y';P ^ nP). (3.22) 

(^) stands for all necessary integrations, while [. . . ] describes all quantum numbers. The amplitude on the 
LHS of Eq. ( 3.22| ), describes all the enhanced diagrams, while the amplitude on the RHS of this equation, 



corresponds to the splitting of one Pomeron to n Pomerons. The precise meaning of this equation, will 
become clear in the next equation. The convenient form of Eq. ( p. 22 ), has been written [20,26,28] in terms 



of the generating functional of Eq. (3.7), and it takes the form 



CO 



1)™ n d n ZP(Y - Y', u p ) d n Z t (Y', u t ) 



N(Y) =J2 7 n - - V : w ' '^ k=i " " ^ '^ k=i- (3-23) 

n=l ^ 

In Eq. ( p.23 ) we denote by Z v and Z l the generating functions that describe projectile and target respec- 
tively. 

Eq. ( |3.23j ) shows that each dipole with rapidity Y' from the target, can interact with any dipole from 
the projectile, (see Fig. ||) with the scattering amplitude 7. The factor 1/n! in Eq. ( |3.23| ) appears due 
to the identity of Pomerons. Eq. ( |3.23| ) is defined in the kinematic region of Eq. ( 3.19[ ), and has a clear 



physical meaning, being the scattering amplitude of two partons (two dipoles) at low energy Yq < 1/Ajp. 
Eq. ( 3.23[ ) gives a natural generalization of Eq. ( 3.18| ), with an obvious physical interpretation, that it is 



the sum of terms, and each of these terms is the product of probabilities to find n-dipoles in the projectile 
and target, multiplied by the scattering amplitude. Since we are discussing the generating functional that 
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Y 



VWWW Pomerom 



O 



A 



Figure 4: An example of enhanced 
diagrams, that contribute to the 
unitarity constraint in the f-channel. 
Wave lines denote the Pomerons. 7 
is the amplitude of the dipole-dipole 
interaction at low energies (at ra- 
pidity Yq w 1/«s )■ The partic- 
ular set of diagrams shown in this 
figure, corresponds to the MPSI ap- 
proach [26,28]. 



satisfies Eq. ( [3.8P as the initial condition, Eq. ( p. 23 ) gives the sum of enhanced diagrams or, in other words, 



at high energy it leads to a new resulting Green's function of the Pomeron. 

The generating functions for the projectile Z p (Y — Y') and for the target Z t (Y r ) in Eq. ( 3,23j ), satisfy 



a very simple equation that describes the parton cascades, in which a parton can only decay into two 
partons. This equation has the form 

- dZ{ ^ U) = -A P u(l - u) 9Z ^ U \ (3.24) 
oy ou 

The above equation has the solution 

u 1 

Z(y,u) = ; r—jr = T , (3.25) 

vy ' ; u + (l-u)e A ^ 1 + -y R e A py y ' 

where u = 1/(1 + 7_r). Eq. ( |3.25| ) satisfies the initial and boundary conditions of Eq. (|3.8| ) and Eq. (|3.9|) . 
Using Eq. ( gjj) and Eq. (gig) it is easy to show that the amplitude is equal to 



N(Y) = 1-Z(u = l/ 1R ) = Y,(-l) n l R e nAlpY . (3.26) 



71=1 



Eq. ( 3.26 ) sums the 'fan' Pomeron diagrams and corresponds to the mean field approximation (MFA) of 



our problem. 

3.4 High energy amplitude 

Using MFA and Eq. ( p\26| ), we can rewrite Eq. Q3.23Q as 



N MPSI {Y) 



n\ 7 W7 R Wjf (3 - 27) 

i - {^ir^w) nMFA ( y - y '^*) n o mfa (y'^ r )}\^ ., 



7^=0 
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Substituting Eq. (|3.25|) into Eq. (|3727| ) , we obtain that [33,34] 

1 \ 1 



N MPSI 



1 



cxp 



r o, 



i 



(3.28) 



T(Y) J T(Y) ' V"' T(Y) 

where T (0,x) is the incomplete gamma function (see formsphl.epsulae 8.350 - 8.359 in Ref. [36]) and 



T(Y) 



7 e 



(3.29) 



In Ref. [34], the solution given by Eq. Q3.28D , is compared with the exact solution (see Fig. 13 of Ref. [34]). 
It turns out to within a 5% accuracy, that the MPSI approximation describes the high energy behaviour 
of the amplitude. At high energy Eq. fl3.28| ) gives N^ PSI (Y) — > 1. 

As has been mentioned, since Eq. (|3.23| ) satisfies the initial condition of Eq. (p.8j) , it describes, as well, 
Eq. ( |3.28| ), the set of enhanced diagrams, and gives the resulting Pomeron Green's function 



G MPSI / y _ y'\ 



N. 



MPSI 
el 



Y - Y') 



(3.30) 



Replacing the bare Pomeron Green's function by Eq. ( |3.30| ), we can use Eq. dig) , Eq. (|2~15|) and Eq. (|1|), 
to calculate elastic and total cross sections. 



3.5 Diffractive production processes 



The conclusion derived from the previous discussion, is that we cannot use the formulae of Eq. ( p. 13 ) and 
Eq. ( |2,14| ), for calculating the single and double diffraction cross sections. Indeed, these two equations 
describe the diffraction production due to the G-W mechanism, while the sum of all enhanced diagrams 
leads to a new source of diffractive production (see Fig. || for examples of such processes). 



Figure 5: Several examples of the 
Pomeron diagrams that lead to a differ- 
ent source of the diffractive dissociation 
that cannot be described in the frame- 
work of the G-W mechanism. Fig. |^-a is 
the simplest diagram that describes the 
process of diffraction in the region of 
large mass Y-Yi = ln(M 2 /s ). Fig. |- 
b and Fig. ||-c give examples of more 
complicated diagrams in the region of- 
sphl.eps large mass. The dashed line 
shows the cut Pomeron, which describes 
the production of hadrons (see Fig. || -a 
which illustrates this point). 





We use the MPSI approximation to obtain the expression for the additional contribution to G-W 
mechanism. The main idea is shown in Fig. |(| As was discussed in Ref. [29] , in order to apply the MPSI 
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approach to diffractive production, we need to consider the generating function of three variables: w,w and 

Vim 

oo 

Z(w,w;v in ;Y), = ^ P(n,n,k\Y) w n w m v\ n . (3.31) 

n=0;m=0;fc=0 

P(n,n, k\Y) is the probability to find n and m Pomerons in the amplitude and conjugated amplitude 
respectively, while k is the number of cut Pomerons (see Ref. [37]). Z (w,w;Vi n ;Y) has been found in 
Ref. [29] for the parton cascade, with a decay of one parton to two partons, 

sphl.eps Z(w,w,v in \Y-Y M = ln(M 2 /s ) = Y m ) = (3.32) 
we -A P Y ^g-AjpY ( w + ^ _ v . n y-A P Y m 



1 + w(e- A iP Y ™ - 1) 1 + w(e- A P Ym - 1) 1 + (w + w - Um)(e" Ajp Ym - 1) ' 



To find the cross section for single diffractive production, we need to calculate the term which is 
proportional to i>j n , and to replace vi n by 2Ajpww. Indeed, this term means that at Y = Y m we have 
only one cut Pomeron, while all other cut Pomerons at this rapidity have decayed to Pomerons without 
cuts. For simplicity we choose Y' = Y — Ym = Y m (see Fig. ^). Replacing Vi n — * 2Ajpww, means that at 
Y = Ym, the last cut Pomeron splits into two Pomerons (see Fig. |J). Therefore, 

N^ FA {w,w;Y-Y M = HM 2 /s )^Y m ) = 2A P «- — (3.33) 

(1 + (w + w) {e A P Ym - l) 2 ) 





Figure 6: The MPSI approximation 
for the cross section of single diffrac- 
tive production of mass (M 2 ,Y — Ym — 
\n(M 2 /so)). The dashed lines shows 
the cut Pomerons. All other notations, 
are as in Fig. H. 
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The general Eq. |[2f) can be rewritten in this case in the form 



N^(Y,Y m = \n(M'/ So )) = (3.34) 

E 7 n+m T^^ N% FA (w?, Y - Y M = HM 2 /s ) Y m ) \ w=1 , a=l 

n! ml a 11 wP o n wP 



n=l:m=l 



d n n mfa ( w t iY -Y m ) d m N MFA (w* , Y - Y m ) 

x lu,t=1 Wrt 



where 



and 



x N MFA (w\Y -Y m = Y M ) N mfa (w\ Y — Y m = Y M ) } \ wP=w t =i!lP= ^ =1 
A 2 e A P (2y-y m ) 

G(L) = L((L-l) 2 -2)+e 1/L (l + 3L)r (l/L) (3.36) 



L (Y, Y m ) = 7 exp (Ajp(Y - Y m )) {exp (AjpY m ) - 1} . (3.37) 

N^ PSI (Y,Y m = ln(M 2 /s )) describes the differential single diffraction cross section for the production of 
mass M. We can calculate the integrated diffraction cross section 

py m (max) 

N^ PSI (Y;M max ,M min ) = / dy m N^ PSI (Y,Y m = \n(M 2 /s )) 

Jy m (min) 

= l e ^Y (B(L(Y,Y m (max))) - B (L (Y, Y m (min)))) , (3.38) 

D 

where Y m {max) = ha.(M^ lax /so) , Y m {min) = ln(M^ in /so) and 

11 eTT (l/L) 
B(L) = 2 + — - ( 3 - 39 ) 

M max and M m i n denote the largest and the smallest masses which are produced in the diffractive process. 

The expression for the integrated cross section for double diffraction can be obtained directly from the 
unitarity constraint of Eq. Q2,4| ), as the diagrams that describe the elastic and single diffraction cross sec- 
tions, do not contribute to the set of Pomeron diagrams, that describe the exact Pomeron Green's function 
(see Fig. which contains examples of the diagrams that contribute to double diffractive production). 

The unitarity constraint is given by 

2N MPSI = N MPSI + N MPSI^ {3AQ) 

where N^ p stands for the inelastic cross section. It was shown that N^ PSI is equal to N MPSI (2T(y)) 
(see Refs. [29,38,39]). Therefore, the integrated double diffraction cross section can be written in the form 

N^ PSI (Y) = 2N MPSI {T(Y)) - N MPSI (2T(F)). (3.41) 
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Figure 7: Several examples of the 
Pomeron diagrams that lead to the dou- 
ble diffractive production. Fig. [7]-a is the 
simplest diagram that describes the pro- 
cess of double diffraction in the regions 
of large mass Y — Y\ = ln(M^/s ) and 
Y2 = ln(Mf/so)- Fig. ||-b contains exam- 
ples of more complicated diagrams in the 
region of large masses. The dashed line in- 
dicates the cut Pomeron which describes 
the production of hadrons (see Fig. ^|-a) . 



3.6 Our approach 

In our approach we combine the G-W mechanism with the exact Pomeron Green's function of Eq. ( |3.30| ). 
First, we replace the bare Pomeron Green's function G(Y — Y') = exp (— Ajp (Y — Y')) in Eq. Q2.7| ) by 
qMPSI (y - Y') of Eq. (|3730[ ), and, obtain 

fyk (s, b) = 9i g k G% PSI (Y - Y') S {b; mi ,m k ; a' P ln(*/*o) = 0) , (3.42) 

with a profile function S (b; rrii, m k ; a'p ln(s/so)) determined by Eq. ( 2.1 0| ) . 

Using Eq. (^T2|) , Eq. (p^) and Eq. (|1|) with U ijk (s, b) from Eq. figl ), we can calculate elastic and 
total cross sections. These formulae provide a correct description of low mass diffractive production due to 
the G-W mechanism. However, to include the Pomeron interactions in processes of diffractive production, 
we need to change Eq. ( |2.13|) and Eq. ( 2.14[ ). We first introduce (see Fig. ||-a) 

0$ (s, M; b) = 9i gl N^ PS (Y = ln(«/*o), Y m = ln(M 2 /s )) S sd (b; m h m k ) , (3.43) 

where the new profile function S sd (b; mi,m k ) is the Fourier transform of 

' 1 ' ' => S sd {b-m^m k ). (3.44) 



l + q 2 /mfY J (2i) 2 l + ( g -+fc)2/ m 2)2 1 + jfe)2/m?) 2 

We found to within an accuracy of around 5%, that the profile function S sd (b; mi,m k ) can be approximated 
by 

2 

S sd (b;mi,m k ) = —^S(b-,m i ,m k ;a ] pln(s/so) = 0) , (3.45) 

12 7T 

where rh\ = 2 \/3 m\ . 

For the calculation of the integrated cross section of single diffractive production we define 

nt d k {s,M max ,M min -b) = (3.46) 
9i9 2 k N^ f p f s (Y = \n(s/s ),Y m (max) = ln(Ml ax /s ),Y m (min) = ln(M 2 in /s )) S sd (6; m h m k ) . 
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For the integrated cross section of the single diffraction channel, we obtain the following expression which 
takes into account both the G-W mechanism, and the enhanced Pomeron diagrams: 



a diff 



M ma ^ M 2 da sd (s,M) 
M 2_ ~W dM 2 



d z b { a 2 (3 2 [a z e 



(a 2 -p 2 )e- 



2 e 



(3.47) 



+ 



{a 2 QfJ 



diff-Sl! 



?2 n diff -fi 2 



+ 2 a 2 2 nf i [ f e- n 'A] i (3.48) 



where ^fjf are given by Eq. ( 3.46j ). The differential cross section for single diffraction in the region 
M > M min is 



M 



2 da sd {s,M) 
dM 2 



J d 2 b {q 2 fifte -01 ' 1 + (3 2 ftfj e - Q2 - 2 +2a 2 2 fige - " 1 - 2 } , (3.49) 



where fi| d fc are given by Eq. flOSp . In Eq. (|3~47|) and Eq. fljUSp we use tt itk of Eq. ( ^42|) . 
For the double diffractive cross section we introduce 

fig (a, 6) = ft 5fc ^il^ 7 (F) S (6; m 4 , m fc ; a'^ ln(s/ So ) = 0) . 
Using this f^l we obtain for the integrated DD cross section 



&dd 



°°dMf 



M 2 



1 Jo 



°dM$ da sd (s,M 1 ,M 2 ) 



M 2 



dM\ dMl 



d l bl a A p A 



n l,2 n 2 

2 e a - + e 2 



+ {a 4 fifje - " 1 ' 1 + (3 4 ni d 2 e- n ^ + 2 a 2 (3 2 nf d 2 e~ ni ^ 



(3.50) 

(3.51) 
(3.52) 



For the integrated single and double diffractive production, the expressions each contains two terms: the 
first is responsible for G-W mechanism for these processes, while the second originates from the large mass 
diffraction of the enhanced Pomeron diagrams. 

Eq. ( fEIlD , Eq. ( gig ) and Eq. &Wj ) with i|fe (s, 6) from Eq. (|3l2|) and Eq. flOap - Eq. (gjg) give 
the full list of formulae in our approach. We wish to emphasize that our approach is based on the pQCD 
input for calculating the enhanced Pomeron diagrams, and on the G-W mechanism for low mass diffractive 
production. 



4. Results of the fit 

4.1 Cross sections and elastic slope 

We have adjusted the parameters of our model which are listed in Table 2, using the formulae of section 
3.6. The fit is based on 55 experimental data points, which includes the p-p and p-p total cross sections, 
integrated elastic cross sections, integrated single and double diffraction cross sections, and the forward 
slope of the elastic cross section in the ISR-Tevatron energy range. The model gives a good reproduction 
of the data, with a x 2 /d.o.f. ~ 1.25. The quality of description of the experimental data is shown in 
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Fig. Fig. 0, Fig. ||, Fig. ^ and Fig. 10. A significant contribution to x 2 /d.o.f. stems from the uncertainty 
for the value of two single diffraction cross sections, and of the total cross section at the Tevatron. The 
X 2 / d.o.f. in Table 2 is calculated neglecting the contribution of the CDF measurement [41] {o~tot = 80 mb. 
at the Tevatron energy). The important advantage of our approach, is that the model provides a very good 
reproduction of the DD data points. In our previous attempt to describe the DD data [6] within a G-W 
approach, it was necessary to assume a non-factorizable contribution for the Pomeron exchange, resulting 
in marginally acceptable results. 
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P 


a'jp 


9i 


92 


rrt\ 


m 2 


0.335 


0.339 


0.012 GeV' 2 


5.82(0.90) GeV- 1 


239.6(37.27) GeV' 1 


1.54 Ge^ 


3.06 GeF 


A K 


7 


a 'm 


9f 


9f 


^0,1 


X 2 /d.o.f. 


-0.60 


0.0242 


0.6 GeV- 2 


13.22 GeV' 1 


367.8 GeV' 1 


4.0 GeV~ 2 


1.0 



Table 2: Fitted parameters for our model which includes G-W mechanism for diffractive production, as well as 
the diffractive processes that stem from enhanced Pomeron diagrams (Pomeron loops). In our approach we have 
extracted a factor of 7 from the product gigk- Therefore, the value of the Pomeron- hadron vertex is equal to 
9i = .9i(from the table) x ^7, and this value is shown in parentheses. 



4.2 t - dependence of the differential elastic cross section 

The behaviour of the elastic cross section in the region of small t, is characterized by B e i and atot, which are 
included in the set of experimental data that we use for our fit. However, we wish to know the scattering 
amplitude at a relatively high value of t > 0.1 GeV 2 , where the simple exponential t behaviour of the input 



elastic amplitude does not describe the data [6]. In Fig. 11 we plot our prediction with the parameters 
of Table 2 for the ^-behaviour of the elastic cross section at the Tevatron energy W = 1800 GeV. We 
reproduce the data quite well, and it is a considerable improvement over the results obtained in [6]. * 

4.3 Mass dependence of the diffractive cross section 

In Fig. |l^, we plot the cross section of single diffraction as a function of 1 — xi = M 2 /s at the Tevatron 
energy. For this cross section, in the region of high mass, we use Eq. (|3.49| ), while for diffraction in the 
region of low mass, we need to make some assumptions regarding the dependence of this cross section on 
M. Following Refs. [7, 42] we assume that the main contribution for the G-W part of single diffractive 
production stems, from the 1FUP1R term, which does not depend on xl- Therefore, the resulting contribution 
has the following form: 

2 da sd (s,M) 

dM 2 dt [ ' 

d(J Low M / M) da mgh M (s M) 

= M 2 8d dM 2 dt ' " (^^RKterm ; G-W contribution) + M 2 sd ^ ' — - (Eq. (gyp ) 

d -B sd \t\ f^High M N% PSI (Y,Y m ;Eq. (pj)) Low M \ 
Sd \ sd NMPSi {Eq .^) +Gsd }■ 

In our parametrization, the scale of the second term is determined by the G-W mechanism. The fact 
that the JFUP1R term is responsible for xl behaviour, is an additional independent input. However, an 
argument for the MJPM term , is that if Ajp = 0, then both IP IP IP and JRJPM lead to a diffractive cross 
section, which is constant as a function of energy. 



*We thank M. Ryskin who pointed out that the model of Ref. [6], gives a minimum at \t\ ~ 0.1 GeV 2 which contradicts 
the experimental data. 
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W=1800 GeV 

do sd /dtd(l-x L ) 



t = - 0.05 GeV" 




1 -x, 



Figure 12: Dependence of single in- 
clusive cross section on 1 — xl — M 2 /s, 
where M is the mass of the diffractively 
produced system. Data are taken from 
Refs. [41,42]. 



5. Predictions for LHC and Cosmic Rays Energies: 6-dependence of the amplitudes 



Fig. 13 shows our prediction for high energy behaviour of the total, elastic and diffractive cross sections as 
well as of the elastic slope. 

At W=14 TeV (LHC energy): a tot = 92.1 mb, a el = 20.9 mb, a sd = 11.8 mb, a dd = 6.08 mb and 
B e i = 20.6 GeV~ 2 . Comparing these results with the prediction of our previous model (atot = 110.5 mb, 
a e i = 25.3 mb, a sd = 11.6 mb, a dd = 4.9 mb and B e i = 20.5 GeV~ 2 at the LHC energy) we see that the 
above comparison 

leads to a few general observations, on which we will elaborate in the Discussion section. 

1) The predictions of our Ref. [6] for atot and a e i at the LHC are considerably higher than the corresponding 
predictions obtained in the present analysis. The difference between the two sets of predictions grow 
monotonically with energy. Note, that the moderate growth of atot and a e \ obtained in the present study 
continues up to energies as high as 10 5 GeV, and probably higher. See Fig. 13 and Table 3. We attribute 
this behaviour to the fact, that in our fit, a'jp is small but not zero. 

2) These features reflect the fact that whereas our previous model has an output compatible with A e jj ~ 
0.08 up to exceedingly high energies, the present model (in which Pomeron enhanced diagrams play a 
significant role), has an output A e ff which decreases monotonically with energy. 

3) The qualitative features of atot and a e i in our two models are also seen in a sd . However, a dd behaves 
differently, slowly approaching a constant value above LHC energies. 

4) Our B sd slope is approximately two times less B e [. At high energy we have only the results of fit to the 
diffraction experimental data given in Ref. [41]. The values of B sd from this fit are B sd = 7.7 ± 0.6 (4.2 ± 
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Figure 13: The energy dependence of cross sections (Fig. |l3|-a) and elastic slope -B e ;(Fig. [l3|-b) for elastic scattering 
and diffractive production. We have plotted 0.25er tot , so as to show all the predictions on the same figure. 



0.5) GeV~ 2 for W = 546 (1800) GeV. Our values (see Fig. |J-b) are B sd = 8.5(9.9) GeV~ 2 at these energies. 
The large difference in the experimental slopes have been discussed in Ref. [41]. 

To investigate the quantitative features of our present model further, we plot in Fig. 14 the b dependence 
of Ai t k(s,b) and a e i(s,b) output amplitudes, at Tevatron and LHC energies. As seen, ^2,2 and A\p reach 
the black bound at relatively low energies, whereas A\ t i is below the black bound, approaching it very 
slowly. 

As noted in Ref. [6], the output cross sections reach the unitarity bound when, and only when, A\i = 
Ai t 2 = ^2,2 = 1- Since, A\^ grows very slowly with energy, we conclude that the very slow approach to 
the unitarity bound observed in Ref. [6], also occurs in the present model. Despite the fact that the two 
models have very different t dependences, we hope that such structure of the amplitude does not depend 
on the model assumptions, but reflects the principle features of the hadron scattering at high energy. We 
shall expand on this issue in the Discussion section. 



6. Survival probability of diffractive Higgs production 

In the following we limit our discussion to the survival probability of Higgs production, in an exclusive 
central diffractive process. Most estimates of the values of survival probability have been made in the 
framework of G-W mechanism, in two channel eikonal models. A general review of such survival probability 
calculations can be found in Ref. [11]. The general formulae for the calculation of the survival probability 
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Figure 14: Impact parameter dependence of A^k and a e i at different energies. 

for diffractive Higgs boson production, have been discussed in Refs. [10-12]. The structure of the survival 

cc 

(I S 2 2ch |) 



probability expression is shown in Fig. |15|.-a. Accordingly, 

N(s 



DfsV 



where, 



N(s) = Jd 2 h d 2 b 2 [Ei, fe < P\i > 2 < P\k > 2 A^s, h) A k H (s, 6a)(l - Af {{s, (bi + b 2 ))) 
D{s) = Jd 2 h d 2 b 2 \Zi k < p\i > 2 < P\k > 2 A i H (s, h) A k H (s, b 2 



(6.1) 

(6.2) 
(6.3) 



< p\i > is equal to (^proton \ ^i) and , therefore, < p\l >= a and < p\2 >= j3. 

A s denotes the soft strong interaction amplitude given by Eq. ( p.7| ),Eq. ( |2.12j ),Eq. ( |2.13| ) and Eq. ( [2.14) ) . 
The form of An(s,b) has been discussed in Refs. [10,11]. In our model we assume an input Gaussian b- 
dependence for the hard amplitudes. 

A? k = A H {s)T? k (b), (6.4) 
where Ah(s) is an s- dependent arbitrary function which does not depend on i, k, and 

2b 2 

^ffc(^) = tt(r h ) 2 e ( ' Ri ' k) ■ The hard vertices and radii Rf k , are constants derived from HERA J /iff elastic 
and inelastic photo and DIS production [40]. 
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a) 



b) 



Figure 15: Survival probability for exclusive central diffractive production of the Higgs boson. Fig. 15 -a shows the 



contribution to the survival probability in the G-W mechanism, while Fig. [lq-b illustrates the origin to the additional 



facte 



8 



enh 



Following Refs. [11,12] we have introduced in the above, two hard 6-profiles 

A%(b) 
<(&) 





exp | 


' b 2 \ 




v 2B*J 


Vp^d 


exp | 


' b 2 \ 




v 2B*J 



(6.5) 
(6.6) 



The values B^=3.6 GeV 2 and B^=\ GeV 2 , have been taken from the experimental ZEUS data on J/\l/ 
production at HERA (see Refs. [11,43]). 

Using Eq. ( 2.12| )-Eq. ( |2.14| ), the integrands of Eq. ( |6.2| ) and Eq. (6.3) are reduced by eliminating 
common s-dependent expressions, 

N{s) = J d 2 b 1 d 2 b 2 [A H (s,b 1 )A H (s,b 2 )(l-As (b = b 1 + b 2 ))] 2 

d 2 b x d 2 b 2 [(1 - a^s^A^ib^ibz) - a sd (s,b) (<(& x ),4f (6 2 ) 

(6.7) 
(6.8) 

Eq. ([O]) does not give a correct estimate for the survival probability, and should be multiplied by a 
factor ((I S 2 nh I)) that, incorporates the possibility for the Higgs boson to be emitted from the enhanced 



+ (6i)<(6 2 )) -ad^A^A^)} 2 , 
D= f d 2 6 1 d 2 fo 2 [^(6 1 )^(6 2 )] 2 



diagrams (see Fig. 15 -b). Therefore, the resulting survival probability can be written as 

(| S 2 |) = (| S 2 enh |) (Eq. fluD) x (| S 2 ch |) (Eq. fj) . 



(6.£ 
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The first attempt to find (| S 2 nh |) was made in Refs. [34,44], where this factor was calculated neglecting 
the fact that the Higgs boson could be produced from the two gluon scattering with a difference in rapidity 
5Yh = ln(Mjj/so) (see for example Refs. [44,45]). The MPSI approach for this case is shown in Fig. [u], 
and it leads to the result: 



(I S 2 enh \)(Y) = ^;N e y PSI (T) = \-T{Y) + e^J (1 + T(Y)) T 1 ' 



T 3 (Y) 



T(Y) 



(6.10) 



where N^f PSI (T) is given Eq. ( 3.2§| ). It was originally suggested to divide this factor by (| 5 
) (Y = SY H ), 



2 

enh 



(i s 2 enh i) on 



(\Sl h \)(Y;Eq.($m) 



(\S 2 enh \)(Y = 5Y H ;Eq. O)' 



(6.11) 



In this paper we take into account 5Yh in a consistent way, which was outlined in Ref. [34], and presented 
in Fig. 17 [16]. At rapidity Y — Y' — hSYjj, one of the partons (Pomerons) will produce a Higgs boson, 
and it should be removed from the cascade evolution. Therefore, those partons which will participate in 
the evolution will be characterized by a new generating function, 



1 



Z [Y -Y' - -8Y H ,u 



(6.12) 



-A F (Y—y'—^6Yu) 9Z MFA (Y - Y' - \SY H ,u) 



1 



' u 



u + (1 



u) e 



-A P {Y~Y'-UY H ) 



2 ■ 



Z (Y — Y' — \SYjj\ u) should be evolved to rapidity Y — Y' using Eq. ( 3.24| ). This evolution results in an 
improved generating function 



Z(Y- Y', u) 



u + (1 - u) e ^ Ap5YH 



(u + (l-u)e A ^ Y - Y )) 



2 ■ 



(6.13) 



Note that for central Higgs boson production \Y = Y' . 



Symmetrically, using Eq. ( |6T^ ) we need to find Z (Y' - \SY Hl u) , and from Eq. ( gTg ), Z (Y', u). The 
result for this function is 



Z(Y',u) 



u + (1 — u) ez 



hA]pSY H 



(u+(l-u) e A iP Y ') 



An,Y'\ 2 



(6.14) 
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Using these generating functions we obtain 
(\S 2 nh (MPSI)\)(Y) = 



where 



v 2 - (-l) n+1 „ d n Z(Y - Y', u p ) . d n Z(Y',u)u t ) 
2^ — ri 7 E^-z \up=i 



n=l 



IV. 



d n vP 



du t 



lu*=l 



= S(T(Y)) - 2e- A ^ Y - SYH)/2 S1(T(Y)) + e- 2AF{Y - 5Y » )/2 S2 (T(Y) ; 
S(T) = 2_ |-T + eT (1 + T) e^To fl 

51(T) = -L |-T(1 + T) + (1 + 2T) T f ± 



S2(T) = -^T (T-l) 2 -2 



i (J !-3/)e^r f i 



r(y) 



7 e 



A^(y-y') _ 1 \ ( e A F r _ A _ 



(6.15) 

(6.16) 
(6.17) 
(6.18) 

(6.19) 




Figure 16: MPSI approach for 
when we neglect SYh = ln(Mjl/so 



02 

D enh 



in the case 



Figure 17: MPSI approach for (| S\ nh |) in the case 
when we take into account SYh — \ti(M h /sq). 



Using Eq. fl6.9| )-Eq, ( |6.18| ), we calculate the survival probability (| S 2 |) for exclusive Higgs production 



in central diffraction. Our results are plotted in Fig. 18. In the following we focus on a detailed discussion 
of our LHC predictions based on the above. Calculating in the framework of a two amplitude eikonal 
model, (| S 2 ch |) for exclusive Higgs production in central diffraction is equal to 2.35%, which is close 
to the value of this part of the survival probability estimated by the Durham group [5,7], and by our 
group [11]. However, the additional factor (| S 2 nh |) = 0.063 as derived from Eq. (6.15), leads to a final 
(| S 2 |) = 0.15%. This result reflects a tendency for the value of the survival probability to be much 
smaller, than when evaluated in the two channel models [6,44]. Note, that the approximate expression for 
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Figure 18: Energy dependence of centrally produced Higgs survival probability. 

(| S 2 nh |) derived from Eq. ( |6.10| ) with our parameters, turns out to be three times smaller. In Ref. [7] 
the Pomeron diagrams were taken into account, but nevertheless, their contribution in the calculation of 
(| S 2 |) was omitted, without any explanation. 



7. Discussion 



In this paper we have presented an approach for soft interactions at high energies based on two ingredients: 

1) The Good- Walker mechanism for elastic and low mass diffraction. 

2) Pomeron enchanced contributions which leads to the exact Pomeron Green's function, which is signifi- 
cantly different from one Pomeron exchange. This component provides the main contribution to high mass 
diffraction. 
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Tevatron 




LHC 




W = W 5 GeV 






GLM 


KMR 


GLM 


KMR 


GLM 


KMR 


a t ot( mb ) 


73.29 


74.0 


92.1 


88.0 


108.0 


98.0 


a d (mb) 


16.3 


16.3 


20.9 


20.1 


24. 


22.9 


a sd (mb) 




9.76 


10.9 


11.8 


13.3 


14.4 


15.7 




Jow M 


8.56 


4.4 


10.52 


5.1 


12.2 


5.7 




high M 


1.2 


6.5 


1.28 


8.2 


2.2 


10.0 


a dd (mb) 


5.36 


7.2 


6.08 


13.4 


6.29 


17.3 


((T e l + &sc 


+ Odd) latest 


0.428 


0.464 


0.421 


0.531 


0.412 


0.57 


S 2ch(%) 


3.2 


2.7 - 4.8 


2.35 


1.2-3.2 


2.0 


0.9 - 2.5 


S enh( % ) 


28.5 


100 


6.3 


100 


3.3 


100 


S 2 (%) 


1.2 


2.7 - 4.8 


0.21 


1.2-3.2 


0.05 


0.9 - 2.5 



Table 3: Comparison of the GLM ( this paper) and KMR [7] models. 



Our enhanced Pomeron formalism, is based on the observation that the soft scattering cross sections 
and slopes, can be reproduced with a'jp « 0.01 GeV~ 2 , rather than a' P = 0.25 GeV~ 2 typical of 
conventional Regge phenomenology. Our result of a very small otjp, implies that the observed shrinkage of 
the forward differential cross sections, traditionally associated with a'jp > 0, can be also reproduced with 
a'jp ~ coupled to strong screening, which produces the desired shrinkage. 

In the following we discuss further the properties of our model, by comparing it with the KMR 
model [7], which is conceptually similar to ours, having the same two mechanisms: G-W and multi- 
Pomeron interactions. In the KMR model a'jp = is an input assumption. In the present GLM model 
a'jp is a fitted parameter, and its value a'jp = 0.012 is an output. As we have discussed the small value 
of a'jp, led us to our key hypothesis that we advocate in this paper: the soft processes are not so soft, but 
stem from short distances, where the QCD coupling is small (ag ~ 0.12 to 0.16). Using this hypothesis we 
built our theoretical approach. This approach is self consistent based on pQCD. This enables us to restrict 
our summation only to triple Pomeron vertices. KMR summation is based on an ad hoc assumption 

r (n Pomerons — > m Pomerons) oc nm\ rn+n ~ 2 . (7.1) 

Further, KMR claim that Eq. ( |7.1| ) leads to their main equation (Eq.(26) of Ref. [7]), which corresponds 
to the parton model. Both claims maybe correct, but they have not been proven in the KMR paper. A 
formal difficulty also noted by the KMR authors (at the beginning of their section 4.3) is that the formulae 
for diffractive production are ad hoc, and are not actually compatible with our Eq. ([O]) and Eq.(26) of 
Ref. [7). 

In spite of the fact, that we do not think that KMR model is able to provide reliable estimates, it is 
interesting to compare our results, since both are based on the same physics. 

1) The introduction of Pomeron induced interactions to the calculation, results in the accumulation of 
Pomeron loops along the initial Pomeron propogator, which lead to a monotonic reduction of the output 
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A e jj with energy. In the KMR model, where additional diagrams, and not only diagrams with triple 
Pomeron interactions have been included, this process occurs more rapidly than in our model. Accordingly, 
they choose a very high input value, Aj n = 0.55. In our fit we have A, n = 0.335. Both in the KMR model 
and in our model the effective shrinkage of the diffraction peak, stems from the Pomeron interactions. The 
difference is that R 2 (s) in the GLM model grows as ln 2 s with a coefficient proportional to a' p . Since in 
evaluating our summations, we have made an approximation in which a'jp = 0, our calculations are only 
trustworthy up to W « 10 5 GeV. The high energy output of GLM and KMR models are presented in 
Table 3. For add the contribution to the diffractive channels coming from the Pomeron enhanced diagrams, 
are larger in KMR than in GLM. Note that add, as calculated by GLM, saturates just above LHC energies, 
while in KMR it continues growing even at energies of W = 10 5 GeV, where it is predicted to be much 
larger than the a s d- 

2) In Table 3 we define c^ ow as the contribution of G-W mechanism, while in the KMR model low and 
high mass diffraction is allocated to mass values M < Mq = 2.5 GeV and M > Mq = 2.5 GeV. 

3) The behaviour of the ratio Rd = i^sld^sd+^dd) conve y S information regarding the onset of unitarity 
constraints at high enough energy. In the G-W model the Pumplin bound Rd < 0.5 is relevant. The 
multi-Pomeron induced contributions are not included in this bound. In the GLM model, Rr> < 0.5 and 
decreases very slowly with energy. It's corresponding R^~ w ~ 0.35, is a constant similar to the output of 
Ref. [6]. In KMR, Rp > 0.5 and grows with energy. The difference between GLM and KMR appears to 
be due to the different summations of the Pomeron induced diagrams, incorporated in the models. 

4) The difference in the values of the survival probabilities calculated in the two channel models are 
significant. The KMR estimate should be reduced by S 2 nh which is ignored in their calculation. Including 
this factor brings the GLM and KMR numbers closer. 

5) The data analysis aimed at determining the three opacities fi^fc, is conceptually different in the GLM 
and KMR models. In principle this information should be obtained utilizing Eq.(2.12)-Eq.(2.18). However, 
the experimental information we have at the UA(4)-Tevatron energies is not sufficient to constrain the 
Pomeron parameters. In the GLM model we, thus, included both the Pomeron and Reggeon trajectories, 
covering also the extensive ISR data. This enabled us to extract the Pomeron (fitted) parmeters. KMR 
have adopted a different strategy, obtaining the Pomeron opacities from a good fit to ^f - in UA(4)-Tevatron 
energy range. This implies a good reproduction of a tot-, ®ei and B e [. From their paper [7], it is not clear 
whether KMR also fit the diffractive channels. This is not the only ambiguity in the KMR presentation of 
their results. Some of their parameters are explicitly presented and fitted, some are presented and "tuned", 
some are assumed and some are implied, but not explicitly presented. 

6) As noted, we attribute our dynamical result to the output of our fit, in which a'jp is small but not 
zero! Both our model and KMR, predict total and elastic cross sections which are significantly smaller 
than in two channel eikonal estimates, which do not include the Pomeron enhanced contributions (see 
Refs. [4-6]). The differences are sufficiently large, so that measurements at LHC and Auger should be able 
to discrimenate between the various approaches. See details below. 

7) The most practical, and perhaps the most interesting result we have obtained, is the small value (about 
0.15%) for the survival probability of central diffractive Higgs production at the LHC. The very small 
value of the final S\\, is due to the smallnes of S 2 ^, multiplied by the small S 2 nh . Our calculation does 
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not include further reductions of S'jj, due to additional short distance processes (see Refs. [44,46]). These 
may further reduce the result we obtained for the survival probability. 

The region of applicability for our formulae is given by Eq. Q3.19 ), in which, Y -C I/7 = 41, and by 
the fact that in our procedure for summing Pomeron diagrams we considered a'p = 0. The first restriction 
leads to a very large kinematic range for rapidity. The seco 

that a'jp ln(s/so) < 0.25/mf in our parameterization, leads to Y < 25 . This region, is marked as 
questionable in Fig. [l3|-a. The set of formulae in section 3.6 do not include the contribution of the 'fan' 
diagram to the elastic amplitude. In our context, it is important that the contribution of such diagrams 
which determines the high mass diffraction, turns out to be small. In the following we list a few experimental 
signatures, which should be measured at the LHC in the near future. These will give us a clue regarding 
the veracity of the models discussed in the paper: 

i) Measurements of atot and <J e i should serve as a critical test for the relevance of Pomeron enhanced 
diagrams. The difference in predictions for models including (excluding) Pomeron enhanced diagrams, 
becomes even more significant at Cosmic Ray energies. The Auger experiment, where we expect results for 
cross sections in the near future, at energies in the 10 5 GeV range, should allow us to discriminate between 
the alternative approaches. 

ii) The Pomeron enhanced contribution to the diffractive channels, as calculated by KMR, is considerably 
larger than our predictions. This is significant for a^d, which acquires a large value in the KMR approach. 

iii) An early estimate of the value of Sfj should be obtained by an LHC measuremnt of the role of hard 
central LRG dijet production in a GJJG configuration, when compared to the pQCD prediction. 

To summarize, we developed an approach which is self consistent, and is based on a perturbative QCD 
input. We hope that a more microscopic approach with roots in QCD saturation, can be built and we 
believe that this paper will contribute to such an effort. 
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